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Abstract 

We analytically and numerically study the analogue of the Parker (magnetic buoyancy) insta- 
bility in a uniformly rotating plasma screw pinch confined in a cylinder. Uniform plasma rotation 
is imposed to create a centrifugal acceleration, which mimics the gravity required for the classical 
Parker instability. The goal of this study is to determine how the Parker instability could be un- 
ambiguously identified in a weakly magnetized, rapidly rotating screw pinch, in which the rotation 
provides an effective gravity and a radially varying azimuthal field is controlled to give conditions 
for which the plasma is magnetically buoyant to inward motion. We show that an axial magnetic 
field is also required to circumvent conventional current driven magnetohydrodynamic (MHD) in- 
stabilities such as the sausage and kink modes that would obscure the Parker instability. These 
conditions can be realized in the Madison Plasma Couette Experiment (MPCX). Simulations are 
performed using the extended MHD code NIMROD for an isothermal compressible plasma model. 
Both linear and nonlinear regimes of the instability are studied, and the results obtained for the 
linear regime are compared with analytical results from a slab geometry. Based on this comparison, 
it is found that in a cylindrical pinch the magnetic buoyancy mechanism dominates at relatively 
large Mach numbers (M > 5), while at low Mach numbers (M < 1) the instability is due to 
the curvature of magnetic field lines. At intermediate values of Mach number (1 < M < 5) the 
Coriolis force has a strong stabilizing effect on the plasma. A possible scenario for experimental 
demonstration of the Parker instability in MPCX is discussed. 
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I. INTRODUCTION 



The Parker (or the magnetic buoyancy) instability arises in a stratified plasma, which is 
partially supported against gravity by a magnetic field. It was originally proposed by Parker 
that the magnetic buoyancy resulting in this system is a mechanism for the formation of 
sunspots [1]. Since then the theory of the magnetic buoyancy instability has been developed 
in numerous papers [2-8] and extended to include various effects, such as pressure anisotropy 
[9], double diffusivity [10-14], rotation [11-13, 15, 16], microturbulence [16], line tying [17] 
and nonuniform gravity [18]. Full 3-D magnetohydrodynamic (MHD) simulations of the 
magnetic buoyancy instability have been performed in the solar context [19-24], including 
some that explore the effects of double diffusivity [25]. At present time this instability is 
accepted in the astrophysical community as the driver of magnetic activity in objects ranging 
from stars to accretion disks and galaxies [26]. In the Sun, the Parker instability may even 
play a role in the global-scale magnetic dynamo [23, 27-29]. 

Despite its important role in understanding the dynamics of the astrophysical magnetic 
fields, no direct experimental investigation of the Parker instability has been performed so 
far. The difficulty is in creating the proper conditions for this instability and, especially, 
in imitating relatively strong gravity in plasma. Nonetheless, the Parker instability could 
be realized in experiments with rapidly rotating plasmas, where the role of effective gravity 
is played by centrifugal acceleration. This possibility of exciting the Parker instability was 
pointed out in theoretical studies devoted to centrifugally confined plasmas [30, 31]. How- 
ever, the instability was not observed in the corresponding experiment. The reason lies, 
perhaps, in the profiles of magnetic field and rotation, which are optimized for plasma con- 
finement and, therefore, for suppressing the instabilities. Thus, the experimental observation 
of the Parker instability is still lacking. 

Our present paper is motivated by the recent results from the Madison Plasma Couette 
Experiment (MPCX) showing a controllable plasma rotation [32]. The MPCX has been 
specifically designed to study phenomena associated with plasma flows [33]. A unique ex- 
perimental setup is implemented in MPCX: a multicusp magnetic field localized near the 
walls of the cylindrical vessel provides plasma confinement and along with the applied elec- 
tric field from the wall electrodes drives the prescribed azimuthal flows (Fig. 1). Different 
types of flows potentially achievable in the MPCX can be used for investigation of such 
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FIG. 1: Madison Plasma Couette Experiment (MPCX): (a) sketch and magnetic field amplitude; 
(b) typical radial profile of cusp magnetic field; (c) electrode configuration near wall. Rings of 
permanent magnets of alternating polarity line the inside of the cylinder with their poles oriented 
normally to the walls. Electrodes are placed between the magnets. The azimuthal velocity at 
the boundary can be adjusted through variation of the E x B drift by changing the voltages 
between electrodes. Helmholtz coils (HC) are used to induce external (mainly axial) magnetic 
field. Removable center core (CC) is shown, but it will not be used in the Parker experiment. 
Figure is adopted from Ref. [35] 

phenomena as magnetic dynamo [34] and magnetorotational instability [35]. 

The goal of this theoretical study is to determine the plasma parameters, flows and 
magnetic fields required for excitation of the Parker instability in a cylindrical geometry and 
to demonstrate the possibility of obtaining such instability in the MPCX. Consideration 
is performed in the frame- work of isothermal compressible magnetohydrodynamics (MHD), 
which appears to be a good approximation for the MPCX plasma. The numerical results 
reported in the paper are obtained using the extended MHD code NIMROD [36], which 
can accurately model both linear and nonlinear plasma dynamics in a specific geometry for 
realistic experimental conditions. The NIMROD code has been used previously to study 
dynamo action [34] and magnetorotational instability [35] in MPCX. 
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In the paper we follow the convention of Ref. [26] and use the term "Parker instability" 
to denote an undular mode of magnetic buoyancy instability, in which the wave vector has 
a component parallel to the magnetic field and, therefore, the perturbed magnetic field 
lines are bent. This mode is opposed to an interchange mode, in which the wave vector 
is perpendicular to the magnetic field and perturbations do not bend the field lines. It 
turns out that in a plasma screw pinch with rigid-body rotation the stability boundaries 
are determined by the modes with non-zero parallel component of the wave vector. These 
field-bending modes are of primary interest in the paper. 

The detailed analysis of the Parker instability was performed by Newcomb [2]. He con- 
sidered a perfectly conducting plasma in a varying with height horizontal magnetic field and 
a uniform vertical gravitational field [similar to configuration shown in Fig. 2(a) for X > 0, 
with gravity acting in the positive X direction G = Goe x , and magnetic field having the 
only component B y (X)]. The equilibrium force balance is 

where P is the plasma pressure and p is the plasma density. According to Newcomb, the 
system is unstable with respect to the Parker instability if for some perturbation £ x the 
energy integral is negative, 



2 J \^{Kl + Kl) 



dX 



+ ^i 2 + ( G °f -^')ie.n«ft<o. (2) 



Here 7 is the adiabatic index, K y and K z are components of the wave vector of the pertur- 
bation. The first two terms in the integral are always positive and represent the stabilizing 
effect of magnetic field lines bending. The instability can arise only when the third term is 
negative. In an isothermal plasma with adiabatic index 7 = 1 and P = pC^, where C s is 
the average sound speed, Eq. (2) is reduced to 

W ~ Stt J U, 2 + Kl 8X + K v B v^ 2CI dX l&l J d V K U ' {6} 

which means that the magnetic pressure must increase in the direction of gravity for the 
system to be unstable. A simple physical interpretation of this is given by Acheson [13]. 
In equilibrium, the magnetic pressure supports more mass against gravity than it would 
be possible in its absence. This extra mass, expressed in the form of magnetic pressure, 
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enters the destabilizing third term in Eq. (3). The situation is somewhat analogous to the 
Rayleigh- Taylor instability, in which the heavy fluid lies on top of a lighter one [37]. 

As seen from Eq. (3), two key elements are necessary for driving the magnetic buoyancy 
instability in plasma: gravity and magnetic field. Below we outline some specific conditions, 
which these elements should satisfy in the experiment. 

(1) Gravity (rotation). Centrifugal acceleration arising in a rotating plasma is the efficient 
way to mimic gravity in the experiment. Obvious disadvantages of such an approach 
are the possible excitation of unwanted instabilities related to the rotation profile, 
such as Rayleigh instability, Kelvin-Helmholtz instability, magnetorotational instabil- 
ity (MRI). A further consequence of rotation is that we might expect the Coriolis 
force to stabilize some motions. If we assume that the rotation is uniform in the axial 
direction, then the most limiting condition on radial profile of angular velocity Q(R) 
comes from the MRI: to have a stable rotation in the presence of axial magnetic field 
one needs to satisfy d(il 2 )/dR > [38, 39]. We consider here only a case of rigid- 
body rotation with Q(R) = Q = const corresponding to effective centrifugal gravity 
G = illRe r . Such plasma rotation has been successfully obtained in recent MPCX 
runs [32]. 

(2) Magnetic field. Only internal (i.e., induced by currents flowing in plasma) magnetic 
fields can be potentially buoyant. This is because magnetic fields contribute to the 
plasma force balance only if the corresponding Lorentz force is not zero. Evidently, 
externally applied fields do not have this property, so they alone cannot induce Parker 
instability. In addition, the magnitude of the field (magnetic pressure) has to be 
increasing in the direction of gravity. One of the feasible ways to create such magnetic 
field in the MPCX is to drive axial current through the plasma, so that azimuthal 
field is induced. However, the curvature of a purely azimuthal magnetic field makes 
plasma susceptible to the so-called "sausage" instability (with azimuthal mode number 
m = 0) [40]. This is obviously unwanted for the Parker instability experiment, as this 
instability may mask the buoyancy instability. To avoid this complication, one can 
add an external axial magnetic field (by driving current in the Helmholtz coils, see 
Fig. 1); this leads to a screw pinch with a helical magnetic field. Such magnetic field 
is also unstable to parasitic kink instabilities [40] , but as we will show, the regions of 
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plasma parameters corresponding to the kink modes and the Parker instability are well 
separated, and these two instabilities can be easily distinguished in the experiment. 
In the following we consider only a case with uniform axial current density (which is 
consistent with assumptions of uniform resistivity profile in isothermal plasma) and 
uniform axial magnetic field, so the radial profiles of field components are B^R) = 
B a R/a and B Z (R) = B z0 = const, where B a and B zQ are constant and a is the radius 
of the pinch (Fig. 2). Thus, the main object of our study is the plasma screw pinch 
with constant pitch and rigid-body rotation. 

We note here that unfavorable curvature of the helical field resulting from cylindrical 
geometry of the screw pinch can lead to development of unwanted kink instabilities. To tell 
whether the instability is due to kink modes or magnetic buoyancy, we also perform the 
stability analysis of analogous plasma configuration in a slab geometry, where the effects 
associated with the field curvature are absent and only magnetic buoyancy can play a role. 
Then we compare and contrast the results obtained in the two geometries. 

The structure of the paper is as follows. In Sec. II, we describe the model used in our 
study. In Sec. Ill, the Parker instability is investigated in a slab geometry and the stability 
boundaries are obtained. In Sec. IV, the linear stability analysis of a plasma screw pinch 
with rigid-body rotation is performed in a geometry of periodic cylinder and the region 
of parameters appropriate for the Parker instability is determined. In Sec. V, we study 
the more realistic case of a bounded cylinder with inclusion of dissipative effects. With 
NIMROD simulations, we explore the nonlinear dynamics of the Parker instability as well. 
In Sec. VI, we summarize and discuss the possible scenarios for experimental demonstration 
of the Parker instability in the MPCX. 

II. MODEL 

We consider the stability problem for two equivalent plasma configurations in two different 
geometries - slab and cylinder (Fig. 2). This is done to distinguish the Parker instability 
in the cylinder from kink instabilities related to the field curvature, as in our slab geometry 
there is no curvature in the imposed magnetic field. The corresponding equilibrium fields 
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FIG. 2: Geometries of the problem: (a) slab; (b) cylinder. 



and gravity are given by 
Slab : 



V eQ = 0, 



X 



B eg = B a — *3 y -r d z qk z , 
a 

R 



e„ + B 7n e x , G = Q X e x ; 



Cylinder : V eq = fi Re v , B eq = B a —e v + B z0 e z , G = 0. 



(4) 
(5) 



Note that the gravity in slab geometry is chosen in such way to be analogous to the centrifugal 
acceleration in cylinder. In both cases the background plasma is generally stratified; the 
exact forms of the density profiles are determined by the force balance equation and will be 
given in Sees. Ill and IV for the respective cases. 

As a frame-work for our study we use isothermal compressible MHD model, which in 
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nondimensional form is 

_ = _ V .(nv), (6) 

n^- = -n(v • V)v - ^Vn + (V x b) x b + ng + u( V 2 v + ^V(V • v) J , (7) 
db 

— = V x (v x b) + r]V 2 h. (8) 

As the motions are themselves isothermal, no energy equation is solved. This corresponds 
to the limit where thermal conduction is much faster than plasma motions, which is the 
appropriate limit for MPCX. Magnetically buoyant motions in this system correspond to 
the limit of doubly-diffusive motions [13], and as such, these isothermal motions are likely 
to be maximally unstable to magnetic buoyancy. In Eqs. (6)- (8) the unit of length is the 
characteristic size of the system a (half- width of the slab layer, radius of the cylinder), 
r, n, v, b and g stand for normalized time, density, velocity, magnetic field and gravity, 
respectively: 

V A+ p V B aG 

T = — t, n = — , v=— , b = — , g = T72", 9 
a po V A B z0 V£ 



where Va = B 0z /^/4iFpo is the Alfven velocity based on the externally applied magnetic 
field B z0 (vertical in slab, axial in cylinder), p is the average mass density. The thermal to 
magnetic pressure ratio (3, dimensionless viscosity v and resistivity r\ are defined as 



8"7rP P c 
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B 2 zQ ' V PqclVa 11 A-KcraVA 
where Po is the average thermal pressure, p is the dynamic viscosity, a is the plasma conduc- 
tivity and c is the speed of light; these parameters are assumed to be constant in time and 
uniform in space. For complete description of equilibrium configurations we also introduce 
Mach number M and pinch parameter 9: 



C s B z q ' 

where Vo = ciQq is the velocity at the cylindrical boundary (peak driving velocity), C s = 
a/Po/po is the average plasma sound speed with adiabatic index 7 = 1 for isothermal model 
and B a is the boundary value of the magnetic field component, which is induced by vertical 
(axial) current in plasma. 

The MPCX plasma parameters can be varied in a wide range allowing experimentalists 
a great flexibility in choosing the regimes of operation (Table I). The dependencies of the 



non-dimensional quantities used in our study on the plasma parameters are presented below: 
thermal to magnetic pressure ratio 

1 Bl Q BUG] ' liUj 

normalized viscosity (inverse fluid Reynolds number based on Alfven velocity) 

" = -4r = °- 88 7 ^w 1 . (11) 

PooVa a [m] B z0 [G] 7V 1/2 [10 18 m" 3 ] A 

normalized resistivity (inverse Lundquist number) 

r 2 7V 1/2 ri0 18 m" 3 l u 1/2 A 

?7= — g— = 0.019 ^ [1U m [ft \ (12) 
^oaV A a[m]5 20 [G]T e 3/2 [eV] 

Mach number 

M- Vo -0 10 V ^M1 (13 ) 

where A is the Coulomb logarithm (typically A ~ 10 — 20). Eqs. (11) and (12) are derived 
from the Braginskii equations for a plasma with singly charged ions in a weak magnetic field 
[41]. The typical values of these quantities in the MPCX are also presented in Table I. 

Note that the cylindrical equilibrium configuration in Eqs. (5) contains neither the mul- 
ticusp magnetic field nor the boundary layers present in the real experiment. In our study 
we do not focus on the details of the plasma driving and simply assume that the equilib- 
rium rigid-body plasma rotation is given a priori. This assumption is consistent with recent 
observations showing such rotation in the bulk of the MPCX plasma [32]. 

Eqs. (6)-(8) must be supplemented with boundary conditions. Here we assume im- 
penetrable, perfectly conducting walls, so the normal components of the velocity and the 
time-varying magnetic field vanish at the boundary T: 

v n \r = 0, b n \ r = 0. (14) 

These conditions are enough in the ideal MHD cases (when rj — v — 0) considered in Sees. 
Ill and IV. In the dissipative MHD case (when 77 7^ and v ^ 0) considered in Sec. V, two 
additional conditions are required, namely, no-slip condition for the tangential velocity and 
absence of the tangential time-varying electric field: 

v t |r = (ve,)t|r, Et|r = ^(V x b) t | r = 0, (15) 
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TABLE I: Parameters of MPCX 



Quantity 


Symbol 




Value 


Unit 


Radius of cylinder 


a 




0.5 


m 


Height of cylinder 


H 






1 


m 


Peak driving velocity 


V 




- 


- 20 


km/s 


Axial magnetic field 


B z o 




0- 


- 100 


G 


Average number density 


N 


10 16 


-10 17 


m -3 


Electron temperature 


T e 




2 - 


- 10 


eV 


Ion temperature 


Ti 




0.2 


: - 2 


eV 


Ion species 




H, 


He, 


Ne, Ar 




Ion mass 




1 


, 4, 


20, 40 


amu 


Thermal/magnetic pressure ratio 




> 


8.8 


x 10~ 5 




Normalized viscosity 


V 


> 


2.1 


x 10- 4 




Normalized resistivity 


V 


> 


1.2 


x 10" 5 




Mach number 


M 




0- 


- 8.5 





where v eq is the equilibrium velocity, and tangential electric field at the boundary is deter- 
mined only by the corresponding component of the current. In the following sections we 
study the stability of the equilibrium configurations [Eqs. (4), (5)] in the frame- work of the 
isothermal MHD model [Eqs. (6)- (8)] with appropriate boundary conditions given by Eqs. 
(14), (15). 



III. SLAB: IDEAL MHD STABILITY 



First, we consider the ideal MHD stability (no dissipation, v = rj = 0) of a stratified 
magnetized plasma in a slab geometry as shown in Fig. 2(a). In this periodic slab geometry, 
there are no Coriolis forces and curvature effects associated with the imposed magnetic field 
that will be present in cylindrical geometries, so the magnetic buoyancy is the only possible 
destabilizing mechanism. The normalized equilibrium fields and gravity in this case are 

v eq = 0, b eg = 9x e y + e z , g = g(x) e x = ^ M 2 x e x . (16) 
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The layer is bounded in the x direction (—1 < x < 1), and periodic in y and z. The x 
component of the force balance results in an equation for the equilibrium density profile 
n eq (x): 

R dry R „ 

M 2 xn eq - 9 2 x. (17) 



2 dx 2 

Taking into account that the average normalized density is 1, we find the solution to Eq. 
(17) in the form: 

, , / 20 2 \ V2Me M2 * 2 / 2 29 2 
nM = { 1 -f3M 2 ) ^erfi(M/^) + fiM*' ^ 
where erfi is the so-called imaginary error function defined via standard error function erf 

as 



2 f ,2 

erfi (z) = —i erf (i z) = —= I e dt. 

J 



o 

The density profiles corresponding to Eq. (18) are shown in Fig. 3(a) for several values of 
Mach number. 

We linearize Eqs. (6-8) near the equilibrium given by Eqs. (16), (18) and introduce the 
plasma displacement vector £(x, y, z)e~ tuJT to describe the perturbations of the physical 
quantities: 

5n = -V • (n^g), 5v = 5b = V x (£ x h eq ). 

The linearized equation for the displacement £ is 

-uo 2 n eq i = F[^] = -^V<5n + (V x 5b) x b eq + (V x b eq ) x 5b + 5n g. (19) 

Due to periodicity of the slab in y and z, we can assume $, depends on y and z as e lk yv+ lk ^ z an d 
consider stability of modes with different wave numbers (k y , k z ) separately. The boundary 
conditions on £ in the x direction follow from Eq. (14): 



lx \x=±l 



= 0. (20) 



Eqs. (19) and (20) constitute an eigenvalue problem. 

In this relatively simple system, without solving the full eigenvalue problem one can 
obtain the stability criterion analytically using the energy principle [42] . It states that the 
static (without flows) ideal MHD system is linearly stable if and only if the potential energy 
of the perturbation, 

W = ~\f £*-F[£]d 3 r, (21) 
11 
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FIG. 3: Equilibrium density profiles for f3 = 8, pinch parameter 9 = 2 and different values of Mach 
number M: (a) in slab [Eq. (18)]; (b) in cylinder [Eq. (28)]. 

is positive for all displacements £ satisfying the boundary conditions in the problem (pro- 
vided that the operator F[£] is self-adjoint). Here displacements are assumed to be complex 
in general and the star denotes the complex conjugation. 

A general expression for the potential energy of the perturbation with wave number 
k = (k y , k z ) in a slab is 

W = \J ((fn eg + b^|V^| 2 + F 2 |£| 2 + ,F((b eg -r)(V-0-(be g -|)(V-r)) 

+ gn eq (UV ■ O + C(V • 0) + ^l^l 2 ) rf3r > ( 22 ) 

where F — k • h eq = k z + k y 6x characterizes the component of the wave vector parallel to the 
equilibrium field. The functional given by Eq. (22) is self adjoint, so the energy principle 
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applies. 

Assuming that F is not identically equal to zero on the interval — 1 < x < 1, we minimize 
W with respect to £ y and £ z . The system is unstable if for some £ x the minimized potential 
energy is negative, i.e., 

+ F 2 \Q 2 - MWl^l 2 ) d 3 r < 0, (23) 

in agreement with Ref. [2] and Eq. (3). Note that the only destabilizing term in Eq. 
(23) is due to joint effect of gravity and plasma current (or magnetic pressure gradient); this 
instability condition does not depend on the density profile explicitly. For modes with k y = 
(which are the analogue of axisymmetric modes with m = in a cylinder) the instability 
condition is simplified to 

6 2 M 2 >k 2 x + k 2 z , (24) 

where k x is the effective wave number in x direction. According to Eq. (24), the system 
becomes unstable when magnetic buoyancy (the term on the left-hand side) overcomes 
magnetic tension (the term on the right-hand side). 

This analytical consideration is confirmed by the results of the numerical solution to 
eigenvalue problem [Eqs. (19), (20)] presented in Fig. 4(1. a-c). As one can see from Fig. 
4(1. a), the behavior of the marginal stability curves are similar for all values of k y in the 
region of large Mach numbers M > 5. In that region, the scaling for the critical (required 
for the instability) value of the pinch parameter is 

9 cr oc M~\ (25) 

which is in agreement with Eq. (24). Thus, we can conclude that such scaling is a universal 
signature of the Parker instability in a slab geometry with equilibrium given by Eqs. (16), 
(18). As follows from Fig. 4(1. b), when 6 > 9 cr the instability develops on the Alfven time 
scale tA = cl/Va [growth rate 7 = Im u is normalized according to Eq. (9)]. From Fig. 4(1. c) 
we see that the critical pinch parameter 6 cr decreases for higher wave numbers k y (shorter 
wave lengths). This result appears only in ideal MHD consideration; as we show below, in 
dissipative MHD system the modes with shorter wave lengths are stabilized. 



Iff F 
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FIG. 4: Results of numerical solution to eigenvalue problems for /3 = 8 in (I) slab with k z = —tt/2, 
(II) periodic cylinder with k z = —tt/2 and (III) bounded cylinder with viscosity v = 0.01 and 
resistivity r\ = 0.01: (a) marginal stability curves on the plane of Mach number-pinch parameter 
(M — 9); (b) growth rate of instability 7 = Im cj as function of pinch parameter 6 for Mach number 
M = 6; (c) critical pinch parameter 9 cr as function of k y (in slab) or m (in cylinder) for different 
Mach numbers M (indicated next to the curves). 



IV. PERIODIC CYLINDER: IDEAL MHD STABILITY 

In this section we study an ideal MHD stability (no dissipation, v = rj = 0) of a rotating 
plasma in a cylindrical screw pinch geometry shown in Fig. 2(b). In this geometry of 
a periodic cylinder, the Coriolis forces and curvature forces associated with the imposed 
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magnetic field are self-consistently included. The normalized equilibrium fields and gravity 
in this case are 

v eg = u re v = \ j^Mr e v , b eq = 9re ip + e z , g = 0. (26) 

Plasma is bounded by a rigid, perfectly conducting cylindrical wall at r = 1. To simplify 
the analysis we assume that the cylinder is periodic in z (axial) directions. The equilibrium 
density is determined from the force balance equation: 

13 dUe<1 = ^M 2 rn eq - 26 2 r. (27) 



2 dr 2 



The solution is 



MV W 2 \ e MV / 2 W 2 
Ueq[r) ~ 2 V [lM 2 )e M2 / 2 -l + /3M 2 ' ( 8) 

where we have taken into account that the average normalized density is 1. Samples of the 
density profiles given by Eq. (28) are shown in Fig. 3(b). 

We linearize Eqs. (6-8) near the equilibrium given by Eqs. (26), (28) and consider all 
perturbations in a reference frame rotating with equilibrium angular velocity co = \J (3/2 M. 
Since the system is periodic in azimuthal and axial directions we take the dependences of 
perturbations on (p and z as e im ^ lk -^ z an d study the stability of modes with different (m, k z ) 
separately. The transition to a rotating reference frame mathematically corresponds to 
changing the eigenfrequency of a mode with given m according to the rule: Co = u> — imu , 
where u and uj are the eigenfrequencies in the rotating and non-rotating reference frames, 
respectively. Introducing the plasma displacement vector £(r, ip, z)e~ WT we then obtain 

Sn = -V • (n^), 5v = Sh = V x x b eq ), 

i.e., the expressions for perturbations of the physical quantities in terms of £ are exactly 
the same as in a static case. In the rotating reference frame, the linearized equation for the 
displacement vector £ is now: 

-u) 2 n eq € + 2n eq u> x5v = F[£] = -^V5n+ (V x 5b) x b + (V x b ) x 5b + 5n u 2 Q r e r , (29) 

where oj> = Wo e z- The boundary conditions on £ follow from Eq. (14): 

£r|r=l = 0. (30) 
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Comparing Eqs. (19) and (29), we see that in a system with plasma rotation there is a new 
effect due to the Coriolis force (the term with 5v), while the centrifugal acceleration plays 
the role of gravity. 

Due to the presence of the Coriolis force in Eq. (29) the energy principle does not give 
a stability criterion in this case. However, as shown in Ref. [43] it still can be applied to 
obtain a sufficient stability condition. Neglecting the Coriolis term, we write the potential 
energy of the perturbation with wave vector k = (m/r, k z ) as 



W 



= \ J ((fn eg + b^|V^| 2 + F 2 ||| 2 + ? F((b eg -r)(V^)-(b eg .|)(V.r)) 
+ f ^ftr(V • C) + C(V • 0) + ^^%| 2 + 2i0Ffe£ - d\ (31) 

where F — k • b eq = k z + m6. Assuming that F ^ 0, we minimize W with respect to £ v and 
£ z . The result is that the system can become unstable if for some £ r 



W min — 2 I l 



+ F 2 |£ r | 2 - 2M 2 ^V|e r | 2 2 V \ir\ 2 j rf 3 r < 0, 

(32) 

where /c 2 = m 2 /r 2 + k 2 z . Note that there are two possible destabilizing mechanisms now: 
one of them is the same as in the slab geometry, and is due to the joint effect of rotation 
and plasma current. This is the Parker instability. The second instability, absent in the slab 
geometry, is solely due to the plasma current, and can exist in the system without rotation. 
These are "sausage" or kink instabilities, depending on the azimuthal wave number m. For 
axisymmetric modes with m — 0, Eq. (32) yields the instability condition in the following 
form: 

2# 2 (M 2 + 2) > k 2 + k 2 z , (33) 

where k r is some effective radial wave number. One can see from Eq. (33), that at large 
Mach numbers M»l the critical value of the pinch parameter scales as 9 cr oc 1/M, which 
is identical to Eq. (25). This suggests that the Parker destabilization mechanism dominates 
in the cylindrical pinch when M ^> 1 if the Coriolis force is ignored. 

To account for the Coriolis force we solve numerically the full eigenvalue problem given 
by Eqs. (29) and (30). The results for modes with different azimuthal mode numbers m 
are presented in Figs. 4(11. a-c) and 5. It appears that the Coriolis force plays a significant 
stabilizing role in the range of Mach numbers 1 < M < 5. This stabilizing effect is amplified 
even more with the increase of (5 (Fig. 5). We also note the existence of unstable "windows" 
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on the plane of parameters (M, 9) for modes with m > 2. Such unstable "windows" have 
been found in a similar screw pinch configuration with rigid plasma rotation in Ref. [44]; 
however, the model used in that paper was not isothermal and equilibrium density was 
assumed constant. Our main finding is that the marginal stability curves for different m 
follow the tension-mediated Parker instability scaling given by Eq. (25) when Mach number 
M > 5 [Fig. 4(11. a)]. At these values of Mach number, the primary destabilizing mechanism 
in the system is due to the Parker instability. 

As one can see from Fig. 4(11. b), for non-axisymmetric modes with m > the growth rate 
of the instability becomes zero at some value of pinch parameter 9 above the critical one, 
$o > @cr- This value is determined by the condition F = k z +m9 = 0, i.e., the corresponding 
mode is a pure interchange mode (perturbations do not bend the magnetic field lines). Such 
mode is a consequence of the assumed periodicity in z; it does not appear in the bounded 
cylinder (Sec. V), where characteristic wave number k z cannot be introduced. Fig. 4(11. c) 
shows that similar to the results in slab geometry, the critical pinch parameter 9 cr decreases 
with increasing azimuthal mode number m when Mach numbers M > 5. 

V. BOUNDED CYLINDER: DISSIPATIVE MHD STABILITY AND NONLINEAR 
DYNAMICS 

In this section we investigate the MHD stability of a rotating plasma screw pinch in the 
more realistic and experimentally relevant geometry of a bounded cylinder [Fig. 2(b)], and 
assuming finite dissipation (here we take v = rj = 0.01). The equilibrium configuration is 
the same as in the ideal case in Sec. IV and given by Eqs. (26) and (28). The stability 
study is performed numerically using the NIMROD code. 

Our first step in this stability study is to solve the eigenvalue problem corresponding to 
linearized Eqs. (6)-(8) with boundary conditions given by Eqs. (14), (15). Using NIMROD 
we solve the initial value problem and determine the growth rate of the fastest eigenmode. 
Since NIMROD uses the Fourier decomposition in periodic ^-direction [36] , and the Fourier 
harmonics are all decoupled in the linearized equations, we are able to determine the growth 
rates for all azimuthal mode numbers m. These results are presented in Fig. 4(111. a-c). 

Fig. 4(111. a) shows the calculated stability boundaries on the plane of Mach number- 
pinch parameter (M — 9). The marginal stability curves for different m follow the tension- 
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FIG. 5: Marginal stability curves in cylindrical system calculated numerically from the eigenvalue 
problem given by Eqs. (29) and (30). Results for different azimuthal modes are shown: (a) m = 0, 
(b) m = 1, (c) m = 2, (d) m = 3. Solid line denotes the case in which the Coriolis force is ignored; 
dashed line is for /3 = 4; dash-dotted line is for /3 = 8. In all cases the axial wave number is 
k z = -it/2. 



mediated Parker instability scaling given by Eq. (25) when Mach number M > 5. In the 
region of small Mach numbers M < 1, the critical pinch 8 cr is practically independent of M. 
Instabilities in this region of parameters are due to the curvature of the magnetic field lines 
and take the form of axisymmetric "sausage" modes or non-axisymmetric kink modes. At 
medium values of Mach number 1 < M < 5, the Coriolis force has a stabilizing effect on both 
types of instabilities. The dependence of the Parker instability growth rate 7 on the pinch 
parameter 9 for Mach number M = 6 and different azimuthal numbers m are presented in 
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Fig. 4(111. b). According to these results, the characteristic time of the instability is the 
Alfvenic time scale tA = cl/Va- 

Fig. 4(111. c) shows the dependence of the critical pinch parameter 9 cr on the azimuthal 
mode number m for several Mach numbers M, which are in the Parker instability region. 
We note that for every M > 5 the critical pinch parameter 8 cr reaches a minimum for 
modes with m = 4 — 5. These modes determine the threshold of the Parker instability. 
This is a consequence of two counteracting effects: the decrease of 9 cr with increasing m 
when dissipation can be neglected [as seen in Figs. 4(1. c) and 4(11. c) for ideal cases], and 
the increase of 6 cr for modes with larger m (shorter wave lengths) due to viscous and ohmic 
dissipation. 

Next, we report the results of 3-D simulations of the nonlinear development of the Parker 
instability. We solve the full system [Eqs. (6)- (8)] with the following parameters: M = 6, 6 = 
1, P — 8, v — 7] — 0.01; at these parameters the system is unstable and the destabilization 
mechanism is dominated by the Parker instability. Fig. 6 demonstrates the time dynamics 
of the kinetic and magnetic energies in the simulations. After the initial linear phase of 
the instability, the energies of the non-axisymmetric parts reach some average level. The 
back reaction of the non-axisymmetric distortions on the initial equilibrium configuration 
leads to its modification. Such modifications in axisymmetric parts of the density and 
magnetic pressure are shown in Fig. 7. As expected, the nonlinear dynamics of the instability 
eliminates its original cause: plasma ( "heavy" fluid) goes down along the effective centrifugal 
gravity, while the magnetic field ("light" fluid) rises up against the gravity. This situation 
in many aspects is similar to the well known Rayleigh- Taylor instability. 

VI. CONCLUSION 

We have performed the stability study of the analogue of the Parker instability (magnetic 
buoyancy) in a rigidly rotating plasma column with constant-pitch magnetic configuration. 
The plasma rotation creates a centrifugal acceleration, which imitates the gravity required 
for the classical Parker instability. To distinguish this instability from instabilities related 
to magnetic field line curvature we have also considered the analogous plasma configuration 
in a slab geometry. 

Applying the energy principle and solving the eigenvalue problem in cylindrical geometry, 
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FIG. 6: Time dynamics of (a) kinetic E^in and (b) magnetic E mag energies of different azimuthal 
modes (labeled) in nonlinear run with Mach number M = 6, pinch parameter 8 = 1, thermal to 
magnetic pressure ratio (3 = 8, viscosity v = 0.01 and resistivity r/ = 0.01. 



we found the marginal stability curves in the plane of Mach number-pitch parameter (M—9). 
It appears that the Parker instability determines the stability thresholds for relatively large 
values of Mach number M > 5; the scaling of the critical (required for instability) pinch 
parameter in this case is 6 cr oc M _1 [Eq. (25)]. The mechanism of the Parker instability has 
some analogy with the Rayleigh- Taylor instability: magnetic field ( "light" fluid) supporting 
plasma ("heavy" fluid) against the gravity is potentially buoyant and tends to rise up. 

At small values of Mach numbers M < 1 the instability in the system is primarily due to 
the magnetic field line curvature ( "sausage" for azimuthal number m = or kink instabilities 
for m/0). 

At medium values 1 < M < 5 the Coriolis force stabilizes both types of instabilities; this 
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FIG. 7: Relative deviations in (r, z) plane of the axisymmetric parts of (c) density and (d) magnetic 
pressure from their respective equilibrium profiles (a) and (b) during the nonlinear phase of the 
Parker instability. Solid black lines in (c) and (d) denote the points where deviation is zero. Arrows 
show the direction of the centrifugal acceleration g cf . Calculations are done with parameters listed 
in caption of Fig. 6. 

effect depends on plasma (3, having a stronger stabilization when is larger. The Coriolis 
stabilization of the kink-like instabilities in the systems with rigid-body rotation is somewhat 
known in the astrophysical community [45] but has not been actively studied in relation with 
experiments. At the same time, this stabilizing Coriolois effect could play an important role 
in modern experimental devices (e.g., tokamaks, reversed field pinches etc.) where plasma 
rotation is commonly observed. 

As follows from our ideal MHD analysis, for a given Mach number M > 5 the onset of 
the Parker instability depends only on the value of normalized gradient of the equilibrium 



21 



magnetic field (pinch parameter 6) but not on its magnitude. This means that the instability 
can be obtained even for infinitely small magnetic field. However, in this case the growth 
rate of the instability approaches zero, because the growth rate scales with the minimized 
potential energy [Eq. (32)] which itself scales as B 2 and thus u 2 oc W min oc B 2 . Therefore, 
to observe the Parker instability in experiment a relatively small magnitude of the equilib- 
rium magnetic field is sufficient; in fact, it is just enough to have a field exceeding the Earth's 
magnetic field. Our analysis has been conducted for isothermal motions in an isothermal 
plasma, as should be appropriate for conditions realized in MPCX. If thermal conductivity 
is decreased, this analysis likely provides a lower bound on the Parker instability. 

Our results suggest that the Parker instability can be achieved in an experiment with 
controllable plasma flows, such as MPCX. Simple estimates show that in order to reach the 
Mach number of M — 6, one may have an argon plasma with ^ — 40, electron temperature 
T e = 4 eV and peak velocity of V = 20 km/s. With an applied axial magnetic field of 
Bo = 4 Gauss, the total axial current required for instability is /total — 1 kA; for a uniform 
distribution of current this will give a pinch parameter of 6 = 1. These parameters may 
soon be achievable in the MPCX. 
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